Zhu-Nakamura theory and the superparabohc level-glancing models 
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We study the applicability of the Zhu-Nakamura theory to a class of time-dependent quantum 
mechanical level-crossing models called superparabohc level-glancing models. The phenomenon of 
a level glancing, being on the borderline between a proper crossing of energy levels and an avoided 
crossing, is also an important special case between the two different approximative expressions 
in the Zhu-Nakamura theory. It is seen that the application of the theory to these models is 
not straightforward. We discuss some possible causes of these difficulties and also compare the 
approximative formulas of Zhu-Nakamura theory to those obtained by the generalization of the 
DDP theory. 



I. INTRODUCTION 

The level-crossing models form a paradigm in the stud 
of quantum dynamics of nonadiabatic transitions [l| 
These models describe quantum systems with coupled 
states for which the corresponding diabatic energy levels, 
i.e, the energies related to the system eigenstates when 
there is no coupling present, depend on some external 
parameter and cross so that the states in question be- 
come degenerate at some parameter values. This type 
of situations appear both in time-independent form, for 
example, in the collisional problems where the nonadi- 
abatic transitions happen effectively at the curve cross- 
ings of potentials that depend on a spatial coordinate 2] , 
as well as in purely time-dependent problems where the 
level crossings are induced by external time-dependent 
fields [3|. The proper frameworks for describing the dy- 
namics in such situations are then the stationary and 
time-dependent Schrodinger equations, respectively. 

The pioneering works on the subject were done by Lan- 
dau, Zener, Stiickelberg and Majorana already in 1932 in 
their studies of the linear crossing problem in which the 
diabatic energies depend linearly on the external coor- 
dinate and the coupling between the diabatic states is 
constant 0-0] • This basic model is nowadays called the 
Landau-Zener (LZ) model. The relevance of this model, 
besides its relative simplicity, comes from the fact that 
the transitions are localized in the vicinity of the cross- 
ings so the above-mentioned behavior of the energy lev- 
els and the coupling often form a good approximation 
to describe the dynamics of many real physical systems. 
This and other level-crossing models have been widely ap- 
plied over the years, for example, to the studies of atomic 
and molecular collisions [l|,[3|, laser-atom interactions Q, 
quantum information processing f*!] and in attempts to 
understand the dynamics of quantum phase transitions 
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One of the recent advances in the field of nonadiabatic 
transitions is the Zhu-Nakamura theory (ZNT) formu- 
lated by Nakamura and co-workers, which the authors 
claim to be a generally applicable and accurate approx- 
imate theory for any level-crossing model and, in this 
sense, to form a complete solution of the problem [l|, [ij] . 
It is based on exact results obtained by Zhu and Naka- 
mura for the linear curve-crossing problem, i.e., the time- 
independent LZ model 13|. By generalizing the coupled 
wave integral method introduced by Hinton [ij], they 
were able to calculate the Stokes constants and thus to 
construct the scattering matrix for the problem. Unfor- 
tunately, the exact analytic expressions for the constants 
are very complicated and therefore not very useful. Be- 
cause of this, they went on, starting from these exact 
results and working initially in the collisional setting, to 
build an approximate theory, i.e. ZNT, by introducing 
several phenomenological corrections to the final analytic 
formulas. 

In the recent years, the purely time-dependent prob- 
lems especially have gained a lot of importance because of 
the progress of experimental methods and laser technolo- 
gies in particular. This has lead to possibilities of control- 
ling the state of a quantum system accurately, for exam- 
ple, by specifically tailored chirped laser pulses [15. 16]. 
Therefore, it is natural that the ZNT has been general- 
ized to handle also these time-dependent problems fl7| . 
This is due to the fact that the time-independent linear 
crossing model (i.e. LZ model) and the time-dependent 
quadratic crossing model (i.e. parabolic model (l8l - l2q ]) 
both have the same analytical structure, namely, they 
both can be reduced to a triconfiuent Heun equation [2l| . 
This allows the same approximate formulas to be used in 
both settings. Despite the many different developments 
that are based on ZNT [22, l2J] and the claim that the 
theory comprises a complete solution to curve-crossing 
problems, there seems to be very few articles in the lit- 
erature studying the basic characteristics of the theory. 
Also, the applicability of the theory must be in princi- 
ple considered for each model independently due to the 
phenomenological character of the theory. 



In this paper, we consider a simple class of time- 
dependent level-crossing models, namely, the super- 
parabolic level-glancing models introduced in [23] • The 
parabolic time-dependencies of the diabatic energies have 
been applied e.g . in studies of laser-induced molecular 
dynamics [1^, [l6| and the parabolic level-glancing model 
was also recently used to study the tunneling between 
different energy bands in the case of merging Dirac cones 
[26| . The level-glancing phenomenon can also be seen as a 
counterpart to the case where the energy matches exactly 
the crossing point energy in the time-independent mod- 
els, a notoriously hard problem to approximate [l|, |25J. 
This particular parameter value remains as an important 
special case also in the ZNT, being the dividing parame- 
ter value between the two approximative formulas in the 
theory. While both of the formulas should be applicable 
to this case in principle, their connection is not entirely 
smooth (see, for example. Fig. 1 in Ref. ^7|) and seems 
to be even less studied. 

We apply ZNT to superparabolic models and discuss 
particular aspects of the theory in this case. We also 
compare it to the other well-established approximative 
theory in the field of nonadiabatic transitions, namely, 
to the Dykhne-Davis-Pechukas theory [28, '29\ studied in 
detail recently in [2J| . The structure of this paper is as 
follows. In the section [ll] we introduce the basic formal- 
ism for both the time-dependent and time-independent 
level-crossing problems and introduce the superparabolic 
level-glancing models. In the section Hill the DDP theory 
and its application to superparabolic models is discussed 
shortly, more details can be found in [2J]. In the section 
IIVI we give a very condensed overview of the ZNT and 
the final recommended approximate formulas in the form 
given in the basic reference [i[. In section |V] we present 
and analyze the results that were obtained by numeri- 
cal calculations and compare these to the approximative 
expressions. Finally, the discussion in Sec. I VII ends the 
presentation. 



II. BASIC MODELS AND FORMALISM 



A. Time-Independent Two-State Processes 



Although we concentrate on studying the time- 
dependent models in this paper, it is important to present 
also some of the main aspects of the time- independent LZ 
model in order to understand ZNT better. Here we give 
only a very minimal overview. For the details, the reader 
is directed to references. 

The model describes a quantum system consisting of 
two states which both experience a potential that is linear 
in the coordinate R and that have a constant coupling 
between them. That is, the system is governed by the 
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time-independent Schrodinger equation 

- [-Fi(i? - Rx) -{E- Ex)] VI = 1^12^2 

(1) 

[-F2{R - Rx) -{E- Ex)] V2 - ■t^i2<^i, 

(2) 

with Fi > 0, V^i2 > and Fi > F2. FiS are the slopes 
of the potentials, ^12 the diabatic coupling and Rx and 
Ex are the crossing point and the energy at the cross- 
ing point, respectively. The case where F1F2 > is 
called Landau-Zener type and F1F2 < is called nona- 
diabatic tunneling type. We consider here only the LZ 
case because only it has a direct counterpart in the time- 
dependent theory. 

By transforming the coupled Eqs. ^ and ^ into a 
momentum representation and redefining the variables 
suitably [i|, Q we can reduce these coupled first-order 
equations to a second-order differential equation for, say, 
the state corresponding to ipi, to a equation of the form 



where 



B"{z) + q{z)B{z)^0, 



f \ 1 -2 |1,'22 7 2n2 

q[z) = la z H — [a z — ) , 



(3) 



(4) 



This differential equation belongs to the class of tricon- 
fiuent Heun equation [2l|. The connection with the pa- 
rameters in the original equations is obtained as 



h^ f (F1-F2) 
2m (21/12) 2 ' 



b'=^lLzM(^E-Ex) 
2FVi2 ^ ' 



(5) 

where it is defined F = \/Fi\F2\ and the independent 
variable is 2 = (2Vi2fc) / F where k is the momentum. 

The important point here is only the physical meaning 
of these reduced parameters. The parameter a? repre- 
sents the effective coupling strength and iP' is the effec- 
tive collision energy. From Eq. ([5]) one can see that 0? is 
always non-negative but that 6^ can be both positive or 
negative, depending on whether the energy E is higher or 
lower than the crossing point energy Ex , respectively. In 
these variables, the celebrated LZ formula for the transi- 
tion probability reads [J| 



PLZ = exp 



4a|6| 



(6) 



Zener obtained this result by reducing the time- 
independent problem described above to purely time- 
dependent one by approximating that the relative nuclear 
motion follows a straight-line trajectory with a constant 
velocity, i.e. R{t) = vt where v'^ = 2{E — Ex) /m,. As 
is well known, this formula works only for energies much 
larger than the crossing point energy Ex- The physical 
reason is, that in a collision process the transition point is 
traversed twice and when E approaches Ex from above, 
the two transitions start to overlap. 



B. Time-Dependent T^vo-State Processes 

Often considered canonical examples for purely time- 
dependent quantum mechanical problems are a spin— ^ 
particle in a time-dependent magnetic field or two elec- 
tronic states of an atom in a chirped or pulsed laser field. 
Another example is the collisional situation, explained in 
the previous subsection, with a heavy enough scattering 
particle so that it can be assigned with a distinct classical 
trajectory. 

In any case, whatever the underlaying physical situa- 
tion is, the coherent dynamics of the two-state system is 
then given by the time-dependent Schrodinger equation 
of the form: 



*>^(v'S-l)l-w- 



(7) 



and if{t) — [ci{t),C2{t)] , where ci(i) and C2{t) are the 
probability amplitudes of the diabatic basis states (fi and 
(^2, respectively. The functions e(t) and V(t) are simi- 
larly called as the diabatic energy levels and the diabatic 
coupling in this time-dependent setting. The crossings 
happen at points of time t^ where e{tc) = 0. 

The Schrodinger equation can be also given in the ba- 
sis of the instantaneous eigenstates of the Hamiltonian 
matrix in ([7]). In this adiabatic basis of the system, its 
energy levels, the adiabatic levels, are given as 



and the adiabatic coupling reads 



7(t) = ± 



V{t)e{t) - e{t)V{t) 



(8) 



(9) 



where the overhead dot stands for time derivative and 
one can fix the sign by fixing the relative sign of the 
basis vectors. In general, the functions e(i) and V{t) 
do not have the same zeros, so the level crossing usually 
appears only as an avoided crossing in the adiabatic basis. 
Moreover, when we have 



\V{t)\<^\e{t)l \t\^±^, 



(10) 



the basis vectors of the two bases coincide, apart from 
the possible swap between the labels, and the initial and 
final probability distributions can be obtained from the 
same expression. In the superparabolic models studied in 
this paper, such a swapping of labels do not occur, and in 
the rest of this paper we take the initial conditions to be 
|c2 (— oo) p = 1 so that both the diabatic and adiabatic 
transition probability, are given by P = |ci (-l-oo) p. 



C. Superparabolic level-glancing models 

The simplest time-dependent model that can take 
into account the double-crossing process in the time- 



independent LZ model as explained previously, is the so- 
called parabolic model for which 



e{t) 



Ae-B 



V{t) = Vo. 



(11) 



The parameters A^ Vo are positive while B can also be 
negative. When this happens, the diabatic energy levels 
do not cross and the transitions are possible only by tun- 
neling while for positive B the energy levels cross twice. 
In the limiting case between the two, namely when _B = 0, 
the levels only touch each other at i = and we call this 
a level-glancing case. 

Obviously, each of these cases are in a direct correspon- 
dence with the different cases in the time- independent LZ 
model. In fact, by making the following correspondencies 
between the variables and the parameters of these time- 
dependent and time-independent models 



t = z, A — a and B ~ b , 



(12) 



and choosing the units in a such way that h — 1 and Vq = 
1/2, we see that differential equation for the probability 
amplitude C2 obtained from Eq. ([7]) with (jlip is reduced 
to an equation that is completely equivalent form as ([3|) . 
The same holds for ci but with the replacement t = —z. 
The parabolic model introduced above was the starting 
point for the ZNT, so it is expected that its approximate 
formulas work well in this case for all the parameter re- 
gions. Here we concentrate on the phenomenon of level 
glancing and set i? = 0. Furthermore, we consider a di- 
rect generalization of the parabolic model, namely, the 
superparabolic models, where the diabatic energies are 
proportional to some even power of time. We can also 
reduce the number of parameters by one, and choose to 
work in units, where fi = \ and the superparabolic level- 
glancing models are defined as [2J| 



e[t) = t 



N 



V{t) ~ a — const, 



(13) 



where N — 2,4, 6, . . . and a is positive. Now, the limit 
a — > is the sudden or diabatic limit while for large a 
the process is adiabatic. As explained at the end of the 
subsection piBp . in both these cases we have P « 0. 

It should be noted that the definition of a is here dif- 
ferent from the one used in [1] or [17] for the parabolic 
model. We denote their parameter as azN and it is ac- 
tually equal to A here. The connection between the two 
parameters is given by the relation 



a^ = 1/(4A). 



(14) 



This difference in notation is of course unfortunate but 
we believe it is better to use the same notation as in 
[2J] for better exposition of the new results and for com- 
parison between the existing ones. Furthermore, this 
should not give rise to additional confusion as the ap- 
proximate formulas of ZNT for the time-dependent case 
have to be transformed from the expressions for the time- 
independent results by replacement of parameters any- 
way. 
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tc satisfying (ITTI) coincide with the zeros of the adiabatic 
levels. As discussed earlier, crossings of the adiabatic lev- 
els are typically avoided ones when i is a real variable, 
so tc is in general complex and usually a branch point of 
the eigenenergy. 

In the case that there are many zero points tc, Eq. 
([T5)) has to be complemented accordingly and there ex- 
ists some rigorous results on the matter [30;, 31]. How- 
ever, as discussed already in the seminal paper of Davis 
and Pechukas 1281 and studied later by Suominen and 
co-workers [3, llj, [S^l, including the contribution of all 
the zero points on the half plane as a coherent sum can 
be very useful in order to obtain a good approximation 
for P when the system parameters are outside the adia- 
batic region. Therefore, we define the generalization of 
the DDP formula as 



FIG. 1. (Color online) Schematics of the time-dependence of 
the energy levels and coupling for the superparabolic models. 
The diabatic levels and the corresponding coupling are drawn 
with solid blue lines while the adiabatic ones are drawn with 
dashed lines in red. 



P, 
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where 



N 

E-p iD(t 

fe=i 



Tfe = 4i lim {t - tl) j{t), 



(18) 



(19) 



III. DYKHNE-DAVIS-PECHUKAS THEORY 

One of the most important results concerning nonadi- 
abatic transitions, and existing prior to ZNT, is given by 
a formula first proposed by Dykhne [29] , which connects 
the structure of the zeros tc of the adiabatic energies 
in the complex plane and the way the adiabatic limit is 
approached by the adiabatic transition probability. This 
was later proved rigorously for a class of two-level Hamil- 
tonians by Davis and Pechukas [281. The Dykhne-Davis- 
Pechukas (DDP) formula is given by 



P 



where 



D{t) 



^-2/m_D(te) 



i£+is)~£^is))ds, 



and tc is defined by the equation 

£+{tc)^£-{tc). 



(15) 



(16) 



(17) 



The method of Davis and Pechukas was to move away 
from the real time axis to the upper half of the complex 
plane and integrate the Schrodinger equation in the adi- 
abatic basis along the level line of D{t) which is defined 
by Im[D{t)] = Im[D{tc)]. The main assumptions al- 
lowing this approach were that there is no crossings for 
real t, that tc is well separated from other zero points 
or possible singularities and that the Hamiltonian is an- 
alytic and single- valued at least in a region of complex t 
plane bounded by the real axis and the level line men- 
tioned above. It is evident from Eq. (|5]) that the points 



and 7(i) is the nonadiabatic coupling defined in Eq. ^. 
In particular, it was found in [24] that this definition is 
necessary to approximate the oscillatory behavior of the 
final transition probability of the superparabolic level- 
glancing models. Indeed, from Eq. (|18p one can be 
see that the existence of multiple zero points tc, k — 
1,2,3,... leads to oscillations in the final state popula- 
tions as the parameters are varied. 



1. Application to superparabolic models 

For the case of superparabolic level-glancing model, 
the generalized DDP formula is obtained and discussed 
in detail in the reference [24J and are simply given below. 
The zero points of the eigenvalues of the Hamiltonian 
defined by Eq. ([13]) are 



^k ^ ^1/N^rni2k-1)/2N^ k^l,2,...,N, (20) 

SO the zero points lie on a circle of radius a^'^ in the 
complex t-plane. The integrals over the adiabatic ener- 
gies are given by 



Dif^c 



rje 



iiT{2k-l)/2N 



where 



and 



ry = 2^A.a(^+i)/^ 



(21) 



(22) 



l^N = 



f ^l^^dy^—B (—,-], (23) 

Jo ^ ^ 27V \2N' 2j ' ^ ^ 



where B{x,y) is the /3 function [3J|. These exphcit ex- 
pressions are also used when discussing ZNT as they form 
the time-dependent version of the phase integrals in that 
context. The factors in Eq. ([T5)) are given by Ffe = (— 1)*^ 



and the points t^ are grouped into pairs with the same 
imaginary part to give the generalized DDP formula in 
the form 



Pddp = 4 



S^ii(-l)''e"'""'"[^ 



(2/0-1)] 



snr 



'7C0S— (2fc-l) 



(24) 



IV. ZHU-NAKAMURA THEORY 
A. Backround and exact results for the LZ case 

Determining the probability of nonadiabatic transi- 
tions can be reduced to calculating the Stokes constants 
as mentioned in the introduction. This way, one obtains 
the connection between the different fundamental asymp- 
totic solutions of the differential equation governing the 
process, that are valid in different regions of the com- 
plex plane [3J|. From these constants, one can construct 
the scattering matrix ('5'(a, 6))^„ (usually called transi- 
tion matrix in time-dependent problems) for the process. 
The transition probability is given in terms of the reduced 
scattering matrix elements (S^)^^ as Pu := | ('S'^)-^2 1^- 
In practice, the Stokes constants are known only for very 
restricted class of differential equations. 

In Ref. 13], Zhu and Nakamura calculated the Stokes 
constants for four different classes of second-order differ- 
ential equations. Among them was the important special 
cases of 



y"{z) + q{z)y{z) = 0, 



(25) 



where q{z) is either a quartic polynomial or a polyno- 
mial where the degree of the highest term is 2n and 
the next highest term is of degree n — 1, where n is 
a positive integer. This means that the solved cases 
contained both the time-independent LZ model and the 
time-dependent parabolic model as well as the super- 
parabolic level- glancing models. 

The important general results of this work were that 
the Stokes constants Ui, i — 1, ...,2{n + l) could be ex- 
pressed in terms of only one of them, say C/i, and that 
this C/i could be expressed as a converging infinite se- 
ries depending on the constants of the polynomial q{z). 
Unfortunately, the expression for Ui is too cumbersome 
and so the results of [l3| are not very transparent to an- 
alytical analysis and therefore of only limited practical 
value (See, e.g., the discussions of their method in [35| 
and i36i])- However, from the exact results one could also 
obtain the general form of the scattering matrix elements 
in terms of the Stokes constant Ui so that for example 
for the LZ type problem (see the subsection lll Ap we have 
the relation 



[Slz) 12 



2iIm(C/i) 



r 



Furthermore, when the transition probability for one pas- 
sage of the crossing point is denoted as 



1 



P = 



1 



\Ui\ 



(27) 



we get the exact result concerning the functional form of 
the overall transition probability as 



Pi2 = 4p(l-p)sin2(V^), 



(28) 



(26) 



where '0 = o.'f'giUi). This form is of course similar to 
the equation for double-passage transition probability de- 
rived by Stiickelberg in the study of atomic collisions [6] . 
In order to overcome the fact that the formulas for 
the Stokes constants are too difficult to be practical, Zhu 
and Nakamura have considered semiclassical approxima- 
tions together with "ad hoc" experimental modifications 
to obtain final formulas for the probability of nonadia- 
batic transitions applicable to general situations. Thus 
the content of ZNT is, in a way, reduced to a set of rela- 
tively simple and compact approximate formulas that are 
obtained by first taking the exact functional form (such 
as Eq. ()28|) ) as the starting point and then making any 
modification to its constitutive parts in order to obtain a 
good approximation. The motivation and guideline be- 
hind the construction of the approximate formulas was 
not only to formulate a general theory that would over- 
come the deficiencies of the previously existing theories 
of nonadiabatic transitions but also to give the final for- 
mulas in a simple and compact form that would avoid 
altogether the complex calculus and contour integration 
and only in terms of the adiabatic potentials. This has 
of course relevance in particular to experimentalists and 
in multilevel situations [13, [S^l ■ 



B. The final recommended formulas of ZNT 

For the direct application of the ZNT, Zhu and Naka- 
mura offer the final formulas to be used in a general level- 
crossing situation, whether time-independent or time- 
dependent. These have appeared in different forms over 
the years, but we take the definitive ones to be those 
given in the appendices of [1|. We consider here first 
the parabolic level-glancing case and then discuss appli- 
cability of the general formulas for the superparabolic 
models. As the parabolic model is the one also studied 



explicitly by Nakamura and co-workers, we expect the 
formulas work well in this first case. Their viewpoint at 
the time was however different as they considered time- 
independent processes and the transition probability for 
a variable collision energy and fixed values of the cou- 
pling strength. On the contrary, we have effectively fixed 
6^ {b = 0) and a^ is related through Eq. (ITil) to the 
independent variable in our case. Also, using that rela- 
tion, it is seen that the parameter range covered here is 
much larger. The superparabolic models are very similar 
in character to the parabolic ones but need to be dealt 
within the context of ZNT by the formulas meant for 
general models, so it is interesting to see how well this 
transition works. 

The final formulas for the time-independent LZ case is 
given for two parameter values separately, when energy 
is higher than the crossing energy (E > Ex) or lower 
than the crossing energy {E < Ex ) ■ As this corresponds 
to either non- negative or non-positive b^ (i.e., the cases 
include equality), we could in principle use either one 
of the final formulas. The formulas for E > Ex are 
somewhat more compact and considered here first. 



1. Double-crossing formulas for the parabolic level- glancing 

model 



The overall probability of nonadiabatic transition is 
now still given by I^E^ but the different terms are given 
by the following expressions: The modification of the 
Landau-Zener formula is given by 



p = exp , , , 

4a \b^ + y/b^ + 0Aa^ + 0.7 

and the phase is given by 



1/2- 



, (29) 



^:^a + (t)s =cr-- + -ln - ]-argT { i- ]--. (30) 

This contains the real and imaginary parts of the phase 
integral 



D{tl) ^(T + i5. 



(31) 



Furthermore, it is advantageous to replace this imaginary 
part 6 in (pO|) by the modification 



5^=\l 



0.8' 



-10" 



(32) 



This modifies the phase ([5(1)) for intermediate and large 
values of a or, equivalently, for small and intermediate 
values of a. Note how the probability of a nonadiabatic 
transition for one passage of the crossing point (I^Hl) differs 
from the original LZ formula ^. It follows, that the ZNT 
could in principle work also in the case 5 = 0. It should 
also be noted that the form of Eqs. ([29l) and ([32)1 are ob- 
tained completely heuristically, i.e., it is not derived from 
anywhere (see, e.g., Ref. [23|). The Zhu-Nakamura the- 
ory actually also offers further approximations to avoid 
the complex integration in pip. However, because the 
phase integral (pij) appears also in the DDP formulas it 
was calculated exactly for any complex crossing point of 
any superparabolic level-glancing model in [2J] and it is 
given by (I2ip . so we do not need those approximations 
here but will discuss about them below along with other 
modifications when considering the superparabolic mod- 
els. Let us now use these results for the parabolic level 
glancing model. 

The transition probability for one passage is now given 
as 



p — exp 



v3/2 



1 



v^ V(0-la-3 -I- 0.7)1/4 



(33) 



For the parabolic level-glancing model, the real and 
imaginary parts of the phase integral are equal as can be 
seen from Eq. (|2ip and we simply have 



(34) 



a = 5 = CO?!'', 
where we have defined the constant 

^^ V?r(i/4^ 

3y2r(3/4) ' 



(35) 



The final transition probability can now be explicitly 
stated 
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Pv 



4e" 
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v^ 
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3/2 
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ca 
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ca3/2 , 
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TT 



(36) 



r 



This can be compared to the DDP result for the parabolic 
glancing model, which now explicitly reads 



P^^p = 4e-2^"'^' sin^ 



v3/2 



(37) 



We can compare the forms of the Eqs. ([57)1 and ([55)1 
and one could think that the argument in the sine func- 
tion of the DDP result is just the first term of the cor- 
responding term in the ZNT result. Furthermore, con- 



sidering that the exact total probabihty is of the form 
P — 4p(l — p) sin^(V') one could guess that the exponen- 
tial in DDP result could be interpreted as p. However, 
although its behaviour is similar to the exponential in 
ZNT result, modifying DDP result this way does not of- 
fer an improvement to the approximation. 



2. Tunneling formulas for the parabolic level-glancing model 

Again, the general formula 

Pi2=4p(l-p)sin2(ar5f/i) (38) 

applies, but now with 

p = ^ (39) 

l + S(ff/7r)e2'^-g2sin2(a) ^ ' 

and the Stokes constant as 



Re(C/i) — cos{a) 



and 



^yB{a/Tr)e'' -gisin^ia) 



(40) 



Im([/i) = sin(cr) 



-2a 



B{(j/TT)e^'' - glsm^{a)coiP{a) 



B{a/TT) 



-,1/2 



2giC0S^(CT) -52 



(41) 



and where the function B(x) is not to be confused with 
P function but is defined as 



Bix) 



2ttx 



2x 



xv^{xy 



(42) 



The functions gi and 172 are again ad hoc modifications 
and are given by 



5i = 1.8(a2)'^"%xp(-<5), 



and 



32 = ^ In (1.2 

■KO 



(43) 



(44) 



and were both actually equal to unity in the original time- 
independent LZ case. The constant b^ does not appear 
in these formulas explicitly but only through the phase 
integrals 5 and a. The final formula for the parabolic 
level-glancing model in this tunneling case is therefore 
obtained just by substituting those terms obtained from 
Eqs. (|M)) and p^ into the above formulas. 



3. Application to superparabolic models 

The above expressions dealt only with quadratic time- 
dependencies, but ZNT can be applied to models with 
general time-dependencies by using these same formulas 
but with replacing the phase inte gral terms a and 5 with 
the phase integral of each model [l^- Also, the diabatic 
parameters (in our case, only the diabatic coupling a) 
should be modified, either by fitting the potential to a 



I 

parabolic one or by using directly the final formulas given 
in LI]: 



a^ = 



and the quantities in the formulas are given as 



and 



Vo^{E2{to)-E,{to))/2, 



[E2{tb) ~ Ei{tt,)][E2{tt) - Ei{tt)] 



[E2{to) ~ Ei{to)Y 



(45) 



(46) 



(47) 



where E2{t) > Ei{t) are the adiabatic potentials, ti, is the 
moment when E2 reaches minimum, tt when Ei reaches 
maximum and to when the difference between the adia- 
batic energies is minimum. The philosophy behind this 
in the ZNT is that the formulas (I45])-(JM1) only refer to 
adiabatic quantities so one could take the experimentally 
measured adiabatic energies Ei and fit them. We can 
of course use the adiabatic model parameters given by 
Eqs. (HI) and (fTS]) . 

In the case of superparabolic level glancing models, 
tb,t,o — 0, so these several points actually coincide and (P 
becomes unity. It is then clear from the above relations 
that the direct application of the final formulas of ZNT 
is not possible because of the several ambiguous "0/0" 
type relations. The diabatic coupling (I46|) . however, gives 
Vq = a as it should. 

The phase integral for the superparabolic models was 
calculated in (PT|) . However, ZNT also offers an approxi- 
mate expression for its calculation in its collection of final 



recommended formulas in [i| as 



N = 2 



'' = -n 



E2{t)dt 



El {t)dt 



(48) 



where 



A 



h-{h + tt)/2 I d? 

dt. 



t + b^ 



(49) 



and it is obvious that the ambiguities arise also here. 

As also the first strategy, namely fitting the super- 
parabolic diabatic energy level to a parabolic one, is 
doomed to fail (for a similar reason why the model can- 
not be linearized at the glancing and dealt with by the 
time-dependent LZ model [i3), and as it is stated in [I4I 
that once a and 6 are obtained the final results are not so 
sensitive to other dependencies of the parameters. Moti- 
vated by this, we take the ZNT formulas of the previous 
subsections with the formula (PT|) to study how well ZNT 
generalizes to the superparabolic models in the next sec- 
tion. 



V. RESULTS 

The resuhs of the DDP theory and the ZNT are 
compared to the results of the numerical simulations in 
Figs. [5] and [31 We take the example cases to be the 
parabolic level-glancing model and the superparabolic 
level-glancing models with A^ = 6 and N = 10. 

As expected, the parabolic model is well approximated 
by the ZNT theory, particularly by the double-crossing 
formula which overlaps with the numerical result. The 
approximate formula derived for the tunneling type tran- 
sitions in ZNT as well as the DDP formula give the right 
characteristic features everywhere and tend to the exact 
result in the adiabatic region. It is also interesting to 
note that the two approximate formulas of the ZNT do 
not coincide. 

When going to the higher values of A^, it is seen that 
the ZNT does not give a good approximation. The tun- 
neling formulas are omitted altogether from the middle 
and bottom plots in Figs. [5] and [3] as they give nonsen- 
sical behavior. For the double-crossing formulas, the ex- 
ponentially decaying part dies of too quickly and there 
is only one peak visible in the linear transition probabil- 
ity plot in Fig. [21 Interestingly however, it is seen from 
the logarithmic plot in Fig. [3l that the oscillations of the 
transition probability seem to have the right frequency, 
so the Stokes phase of the LZ model with the modifica- 
tions of the ZNT seems to work well and the problem 
is in the expression for the probability of one transition, 
i.e., in Eq. (|^ . 

The DDP formula on the other hand generalizes easily 
for the higher superparabolic models (i.e. higher values 




FIG. 2. (Color online) The DDP approximation is the blue 
line, ZNT (diabatic crossing case) is the red and the numerical 
result is the black line. The green line in the first plot is the 
ZNT formula for tunneling case. It does not generalize well 
at all for higher values of A'' and is omitted. 



of A) and apart from the small phase shift is a very good 
approximation in the whole parameter region. 



VI. CONCLUSIONS 

We have reviewed the basic aspects of ZNT and studied 
their application to the superparabolic models. Although 
ZNT has been studied in the context of the parabolic 
model before, we considered this case from a different 
viewpoint than done in the original works of Nakamura, 
Zhu and co-workers. Our treatment was explicitly time- 
dependent and we studied the transition probability as 
a function of the diabatic coupling, instead of the usual 
treatment of the equivalent coUisional model where the 
independent variable is the collision energy and the cou- 



N = 2 




FIG. 3. (Color online) This is similar to the previous figure 
but with logarithmic ordinate. 



this region was of course one of the major motivating 
factors behind the whole Zhu-Nakamura approach, as it 
was not properly dealt with by the previously-existing 
theories. 

The level-glancing case resembles in many ways the 
double-crossing one, as there are oscillations in the final 
populations and the adiabatic coupling has two distinct 
peaks, for example. Indeed, the ZNT approximation de- 
rived for the double-crossing case is very accurate in the 
parabolic level-glancing model. As the above-mentioned 
behavior is common to all of the superparabolic models 
it is somewhat surprising that the ZNT at its present 
does not work well for the models with higher values of 
N but that at least the expression for p (Eq. ([29| ) should 
be modified further. It also seems clear, that the com- 
pact formulation of the ZNT for general models, namely, 
the parameter-fitting formulas (1451) - (|47)) , do not seem to 
be straightforwardly applicable to adiabatic energies that 
are symmetric in time, i.e., when to, tb and U coincide. 

At the same time one can see that the generalization 
of the DDP theory, namely, the full summation formula 
(jlSp . can be directly applied in a systematic fashion and 
it does not lead to difficulties in the conceptual level. 
Of course it also lacks of a sufficient general mathemati- 
cal proof that would include superparabolic models and 
more work remains to be done in this direction. On the 
other hand, it may be more difficult to handle multilevel 
problems within DDP theory [37|- Furthermore, DDP 
theory is not so directly related quantities measured in 
experiments as ZNT due to the fact that it relies on ana- 
lytic continuation of the adiabatic energies. Therefore, it 
has been our point to show the usefulness of the both of 
these approaches and especially to highlight the need for 
more general formulation of the ZNT and also the need 
for more careful instructions to its application. 
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